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Abstract. Grids of models of massive stars (M > 20 Mq) with rotation are computed for metallicities Z ranging 
from that of the Small Magellanic Cloud (SMC) to that of the Galactic Centre. The hydrostatic effects of rotation, 
the rotational mixing and the enhancements of the mass loss rates by rotation are included. The evolution of the 
surface rotational velocities of the most massive 0-stars mainly depends on the mass loss rates and thus on 
the initial Z value. The minimum initial mass for a star for entering the Wolf-Rayet (WR) phase is lowered by 
rotation. For all metallicities, rotating stars enter the WR phase at an earlier stage of evolution and the WR 
lifetimes are increased, mainly as a result of the increased duration of the eWNL phase. Models of WR stars 
predict in general rather low rotation velocities (< 50 km s~^) with a few possible exceptions, particularly at 
metallicities lower than solar where WR star models have in general faster rotation and more chance to reach 
the break-up limit. The properties of the WR populations as predicted by the rotating models are in general 
in much better agreement with the observations in nearby galaxies. Some possible remaining difficulties in these 
comparisons are mentioned. The evolution of the chemical abundances is largely influenced by rotation in all 
phases from the MS phase to the WN and WC phases. We also show that the interval of initial masses going 
through the LBV stage is changing with Z and fi. 

The observed variation with metallicity of the fractions of type Ib/Ic supernovae with respect to type II supernovae 
as found by Prantzos & Boissier 12003t is very well reproduced by the rotating models, while non-rotating models 
predict much too low ratios. This indicates that the minimum initial masses of single stars going through a WR 
phase are consistently predicted. At Z — 0.040, stars with initial masses above 50 AIq reach a final mass at the 
time of supernova explosion between 5 and 7.5 Mq, while a.t Z = 0.004, like in the SMC, the flnal masses of stars 
are in the range of 17 - 29 Mq. On the whole, rotation appears to be an essential parameter even for the WR 
properties. Detailed tables describing the evolutionary tracks are available on the web. 
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1. Introduction 

Wolf-Rayet stars are considered to be bare stellar cores 
whose original H-rich envelopes have been removed either 
by strong stellar winds or by mass transfer through Roche 
Lobe Overflow in close binary systems (Contill976ll Chiosi 
& Maeder [THHi; Abbott & Conti il987f) . Their associa- 
tions with young star forming regions implies that their 
progenitors must be massive stars (see e.g. the recent re- 
view by Massey 120031 and references therein) . WR stars 
have a deep impact on their surroundings thanks to their 
high luminosity and their strong stellar winds. Their broad 
emission lines can be detected in the integrated spectrum 
of remote galaxies (Kunth & Sargent 119811 Schaerer et 
al. I1999|l enabling us to study star formation and evolu- 
tion in very different environments, from metal poor blue 



compact dwarf galaxies to the vicinity of AGN (see e.g. 
Lfpari et al. 12003(1 . They contribute to the enrichment of 
the interstellar medium by newly synthesized elements. In 
particular their winds at high metallicity may be heavily 
loaded with carbon (Maeder 1 1992|l . Their winds may also 
eject significant amounts of ^^Al (see e.g. Vuissoz et al. 
201) l.i , responsible for the diffuse emission at 1.8 MeV ob- 
served in the plane of our Galaxy (Prantzos & Diehl ll9t?Sjl , 
in ^^F (Meynet & Arnould|20001 whose origin still remains 
largely unknown (Cunha et al. 2003 ) , and in s-process el- 
ements (see e.g Arnould et al.|l997). The WC star winds 
are also rich in ^^Ne, which explains the high ^^Ne/^°Ne 
isotopic ratio observed in the galactic cosmic ray source 
material (see e.g. Meynet et al. 200Q. WR stars are also 
the progenitors of type Ib/Ic supernovae (see the review 
by Hamuv .20051 . Recently the spectrum of such a super- 
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nova was observed in the optical transient of a 7-ray burst 
(see e.g. Hjorth et al. 12003(1 . confirming the suspected Hnk 
between these stars and the long 7-ray bursts (Woosley 
I1993|> . For all these reasons Wolf-Rayet stars appear as 
objects worthwhile to be well understood. 

In a previous paper (Meynet & Maeder I2003L paper 
X) , we discussed the consequences of rotation on the prop- 
erties of WR stars at solar metallicity. One of the main 
conclusions is that the theoretical predictions for the num- 
ber ratios of Wolf-Rayet to 0-type stars, for the ratio of 
WN to WC stars and for the fraction of WR stars in the 
transition WN/WC phase, are in good agreement with the 
observations when the effects of rotation are accounted for 
in stellar models. In contrast, the models with present- 
day mass loss rates and no rotation do not succeed in 
reproducing the observed values. The main purpose of the 
present paper is to explore the case of metallicities lower 
and higher than solar and to see if the above conclusions 
still hold. 

Let us recall that interesting questions arise from the 
observations of WR stars both at low and high metallic- 
ity. At low metallicity, it has generally been thought that 
WR stars might preferentially be formed by mass trans- 
fer through Roche Lobe Overflow in close binary systems. 
For instance, it was thought that the majority, if not all 
the WR stars in the SMC should be born thanks to the 
mass transfer mechanism. The main reason is that, at low 
metallicity, the mass loss rates are much lower than at 
higher metallicity, thus making the ejection of the H-rich 
envelope by stellar winds more difficult. However, this idea 
was recently challenged by the works of Foellmi et al. 
H2()()3al I2008b|l . They looked for periodic radial velocity 
variability in all the WR stars in the Small Magellanic 
Cloud and in two thirds of the WR stars in the Large 
Magellanic Cloud. They found that the percentage of bi- 
naries among the WR stars is of the order of 40% for the 
SMC and 30% for the LMC, thus comparable or even be- 
low the percentage of binaries among the WR stars in our 
Galaxy. This means that, in the SMC, at most 40% of 
the WR stars could originate from mass transfer through 
Roche Lobe Overflow (RLOF) in a close binary system. 
The real fraction is likely lower since RLOF has not nec- 
essarily occurred in all these systems. Therefore even in 
the SMC, a large fraction of the WR stars likely originate 
via the single star scenario. 

How is it possible ? Docs this mean that the mass loss 
rates are larger than usually found or that another pro- 
cess is at work which favours the evolution of massive stars 
into the WR phase ? In an earlier work (Maeder & Meynet 
we explored the first hypothesis, namely the effects 
of enhanced mass loss rates. We multiplied by a factor of 
two the mass loss rates given by de Jager et al. 1(198811 
during the 0-type star phase and we adopted an average 
mass loss rate of 8T0~^ Mq yr~^ during the WNL phase, 
i.e. twice the average mass loss rate given by Abbott & 
Conti 1(1987(1 . Doing so, we obtained a good agreement 
with the observed WR populations for the metallicities 
between that of the SMC and up to twice the solar metal- 



licity. However, nowadays clear evidence of clumping in 
stellar winds of massive stars (Nugis et al. 11998(1 has been 
found and the new estimates of the mass loss rates (Nugis 
& LamersEOHni are reduced by a factor 2 to 3 with re- 
spect to the enhanced mass loss rates described above. 
Therefore the enhanced mass loss rate hypothesis is ruled 
out and another process must be at work. In paper X, we 
showed that rotation might well be this process. Indeed, 
in massive star models, rotation favours the evolution into 
the WR phase in two ways, first by allowing the star to 
enter into the WR phase at an earlier stage, thus making 
the WR lifetime longer and secondly by allowing smaller 
initial mass stars to go through a WR phase. These effects, 
as we shall see below, are also present at low metallicity 
and allow us to reproduce the observed WR populations 
in the LMC and SMC without assuming that a large frac- 
tion of WR stars owe their existence to mass transfer in 
close binary systems. 

We shall also study the effects of rotation at higher 
metallicity than solar, more precisely at twice the solar 
metallicity, a value often quoted for the galactic centre 
(although some authors quote for this region a value of 
the metallicity similar to that of the solar neighbourhood, 
see Carr et al. 119991 Najarro 2003). This region of the 
Milky Way is very rich in massive stars (according to Figer 
ct al. 120021 the Arches cluster contains about 5% of all 
known WR stars in the Galaxy) and it is thus interesting 
to derive the properties of the WR stars predicted by the 
present rotating models. Moreover, populations of Wolf- 
Rayet stars at a higher metallicity may also be of interest 
for studying the stellar population in the vicinity of AGNs 
(see e.g. the review by Heckman 11999(1 . 

Sect. 2 briefly summarizes the physics of the models. 
Evolution of the surface rotational velocities is discussed 
in Sect. 3. The evolutionary tracks are presented in Sect. 

4. The effects of rotation on WR star formation and WR 
lifetimes at different metallicities are discussed in Sect. 

5. Comparisons with the observed WR populations are 
performed in Sect. 6. Finally the predicted surface abun- 
dances of the present WR stellar models are discussed in 
Sect. 7. 

2. Physics of the modeis 

Most of the physical ingredients of the present models are 
the same as in the solar metallicity models of Meynet & 
Maeder ((2003| paper X). They differ in only two points: 

— The initial compositions are adapted for the different 
metallicities considered here. For a given metallicity Z 
(in mass fraction), the initial helium mass fraction Y 
is given by the relation Y = Yp + AY/AZ ■ Z, where 
Yp is the primordial helium abundance and AY/AZ 
the slope of the helium-to-metal enrichment law. We 
use the same values as in Maeder & Meynet flUUT^ 
i.e. Yp = 0.23 and AY/AZ = 2.5. For the metallic- 
ities Z = 0.004, 0.008 and 0.040 considered in this 
work, we thus have X = 0.757, 0.744, 0.640 and Y 
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Age [10^ years] Age [10^ years] Age [10* years] 

Fig. 1. Evolution of the equatorial velocities {upper panels) and of the ratio ^ of the angular velocity to the critical 
angular velocity {lower panels) at the surface of star models of different initial masses and metallicities. Values of the 
initial masses are indicated in the upper left and right panels. The initial masses corresponding to the tracks plotted 
in the middle panels are the same as those shown in the left panels. The initial rotational velocity of all the models is 
300 km s~^. Empty circles are placed at the stage when the star enters into the Wolf-Rayet phase. 



= 0.239, 0.248, 0.320 respectively For the heavy el- 
ements we adopt the same mixture as the one used 
to compute the opacity tables for solar composition. 
In that respect the present grid differs also from our 
previous non-solar metallicity rotating models (papers 
VII and VIII). In these models, we adopted the "en- 
hanced alpha elements" opacity tables of Iglesias & 
Rogers p996f) instead of the "solar composition" ones. 
As the relative abundances of a-nuclei vary accord- 
ing to the galactic history, the present choice may be 
better. 

— Secondly the wind anisotropics induced by rotation 
were neglected. This last choice appears justified in 
view of the results obtained in paper X. Indeed for the 
initial velocities considered (wini = 300 km s~^), the 
effects of the wind anisotropics have been shown to be 
very small. Let us however emphasize that this is not 
true for higher initial velocities fMaeder 12002(1 . 



Since mass loss rates are a key ingredient of the models 
in the mass range considered here, let us recall the pre- 
scriptions used. The changes of the mass loss rates M with 
rotation are taken into account as explained in Maeder & 
Meynet H2000a|l . As reference mass loss rates we adopt the 
mass loss rates of Vink et al. f2000| I2001|l who take ac- 
count of the occurrence of bi-stability limits which change 
the wind properties and mass loss rates. For the domain 
not covered by these authors we use the empirical law 
devised by de Jager et al. Ijl988fl . Note that this empir- 
ical law, which presents a discontinuity in the mass flux 
near the Humphreys-Davidson limit, implicitly accounts 
for the mass loss rates of LBV stars. For the non-rotating 
models, since the empirical values for the mass loss rates 
are based on stars covering the whole range of rotational 
velocities, we must apply a reduction factor to the empir- 
ical rates to make them correspond to the non-rotating 
case. The same reduction factor was used as in Maeder & 
Meynet (jjSOOlJ. During the Wolf-Rayet phase we use the 
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mass loss rates by Nugis & Lamers H2UL)U|I . These mass loss 
rates, which account for the clumping effects in the winds, 
are smaller by a factor 2-3 than the mass loss rates used 
in our previous non-rotating "enhanced mass loss rate" 
stehar grids (Meynet et al. ll994|l . 

During the non-WR phases of the present models, we 
assumed that the mass loss rates depend on the initial 
metallicity as M{Z) = {Z/ZqY''^M{Zq) (Kudritzki & 
PulsEOnni Vink et al. 12(1171)1 . For models at Z = 0.040, 
we also compute a series of models with metallicity de- 
pendent mass loss rates during the WR phase. According 
to Crowther et al. I)2002|l mass loss rates during the WR 
phase may show the same metallicity dependence as the 
winds of 0-type stars, i.e. scale with ~ {Z/ZqY/"^. 

The effective temperature of Wolf-Rayet stars is a del- 
icate problem, since the winds may have a non-negligible 
optical thickness. Here we adopt a simple correction 
scheme to take account of this effect (e.g. Langer(T989): 
the effective radius ReS at the optical thickness r = 2/3 
is related to the classical photospheric radius R by the 
relation 

_ mm 

^eff — ri+ — , 

where K is the opacity and the other symbols have their 
usual meaning. More details on the procedure to estimate 

T and v^o are given in Schallcr et al. H1992|l . The effec- 
tive temperature at r = 2/3 is then obtained by the usual 
relation L = A-kR^^cjT^^. Such a correction has been ap- 
plied in the WR stages and only there. 

A moderate overshooting is included in the present 
rotating and non-rotating models. The radius of the con- 
vective cores are increased with respect to their values ob- 
tained by the Schwarzschild criterion by a quantity equal 
to 0.1 Hp, where Hp is the pressure scale height estimated 
at the Schwarzschild boundary. The effect of rotation on 
the transport of the chemical species and of the angular 
momentum are included as in our papers VH and VIII. 

As initial rotation, we have considered a value equal 
to 300 km s"i on the ZAMS for aU the initial masses and 
metallicities considered. At solar metallicity, this initial 
value produces time-averaged equatorial velocities on the 
MS well in the observed range, i.e. between 200 and 250 
km s~^. At low metallicities this initial rotational velocity 
corresponds also to mean values between 200 and 250 km 
s^^ on the MS phase, while at twice the solar metallicity, 
the mean velocity is lower, between 160 and 230 km s~^ 
(see Table nj. Presently we do not know the distribution 
of the rotational velocities at these non-solar metallicities 
and thus we do not know if the adopted initial velocity cor- 
responds to the average observed values. It may be that 
at lower metallicities the initial velocity distribution con- 
tains a larger number of high initial velocities (Maeder et 
al. [1999) , in which case the effects of rotation described 
below would be underestimated at low metallicity. 

All the models were computed up to the end of the 
helium-burning phase. Their further evolution in the ad- 



vanced stages will be presented in a forthcoming paper 
(Hirschi et al. in press). In order to facilitate future de- 
tailed comparisons, we provide electronic tables describ- 
ing the evolutionary tracks presented in this paper^. For 
each initial mass, we extracted 350 points describing the 
whole sequence from the ZAMS to the end of the core He- 
burning phase. At each time step we provide the following 
quantities: age in years, actual mass in solar masses, log 
L/Lq, log Toff (which takes account of the optical thick- 
ness of the winds for WR stars, see above), the surface 
abundances in mass fraction of H, He, ^^C, ^^C, ^"'N, ^^O, 
"O, 1^0, ^°Ne, ^^Ne, the fraction of the mass of the star 
occupied by the convective core, the uncorrected effective 
temperature, i.e. which does not take account of the op- 
tical thickness of the winds (for WR stars only) , the mass 
loss rate in solar masses per year, log pc, the central den- 
sity, log Tc, the central temperature, the abundances at 
the centre of the star in mass fraction of the same ele- 
ments as above, the ratio of the polar radius to the equa- 
torial radius, the surface equatorial velocity in km s~^, 
and the ratio of the surface angular velocity to the break- 
up velocity. For the nuclear burning phases, models with 
the same number in the tables always correspond to the 
same evolutionary stage, i.e. have the same mass fraction 
of hydrogen or helium at the centre: point 1 corresponds 
to the ZAMS stage, point 100 to the end of the core H- 
burning phase, point 201 and 350 to the beginning and 
end respectively of the core He-burning stage. 

To conclude this section, let us emphasize that the 
amplitudes of the effects studied in this paper are partic- 
ularly sensitive to the mass loss rates, the extension of the 
convective cores and the initial rotation. Does this mean 
that the present results are very model-dependent ? We 
do not think so for the following reasons: firstly we are not 
free to change the above parameters beyond certain lim- 
its, for instance the initial rotation was chosen in order to 
reproduce the averaged observed velocity in the MS band, 
mass loss comes in part from the radiation wind theory, 
in part from empirically based relations, the overshoot is 
constrained from comparisons between the observed and 
computed MS width; secondly the present physical ingre- 
dients already allowed us to reproduce many observed fea- 
tures such as the observed surface enrichments at different 
metallicities as well as the the blue to red supergiant ratio 
in the SMC. These good agreements give some support 
to the chosen set of physical ingredients; thirdly we think 
that, in this first extensive study of the effects of rotation 
on the WR formation at different metallicity, it is nec- 
essary to focus on the effects of rotation alone, all other 
physical ingredients being kept the same. In the future, 
improvements of the physical ingredients of the models 
will change some of the results but likely the effects of ro- 



The tables can be found at the web address: 
http://ww'w.umge.ch/sciences/astro/an/ choose on the left 
"RESEARCH GROUPS", then under the title "EVOLUTION 
STELLAIRE", choose "ETOILES MASSIVES", then "Tables 
des resultats des recents modeles avec rotation" . 
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tation discussed here will qualitatively remain the same. 
In that respect it is interesting to mention that in Meynet 
H2(JU(J|I we deduced from a previous grid of Wolf-Rayet 
stellar models at solar metallicity (paper V) the theoret- 
ical value for the number ratio of WR to 0-type stars. 
The results obtained (WR/0 = 0.026 for Vi^i = km s"^ 
and 0.072 for Vini = 300 km s~^) are very similar to those 
obtained from the present grid (WR/0 = 0.02 for Vini — 
km s~^ and 0.07 for Vini = 300 km s~^, see Table EJ, de- 
spite the fact that the two sets of stellar models differ 
in many respects. Among other things, they differ by the 
prescription for the mass loss rates, the amount of over- 
shooting and the form of the shear diffusion coefficient. In 
our view, the nearly identical results obtained in the two 
cases supports the idea that the results concerning the 
WR populations are not very sensitive to small changes of 
these parameters. 

3. Evolution of the rotational velocities at various 
metallicities 

As discussed in paper X, the evolution of the rotational 
velocities at the stellar surface depends mainly on two fac- 
tors, the internal coupling and the mass loss. The internal 
coupling is achieved in radiative zones by meridional circu- 
lation and by shear diffusion. During the Main Sequence 
phase, the internal coupling transports angular momen- 
tum from the contracting core to the expanding envelope 
and thus maintains the surface angular velocity above the 
value it would have in case of no-coupling i.e. of local 
conservation of the angular momentum. Mass loss on the 
other hand removes the angular momentum contained in 
the ejected outer surface layers. Thus, the two effects act in 
opposite directions for the evolution of the surface veloc- 
ities. The stronger the coupling mechanisms, the higher 
the values of the surface velocity at a given stage. The 
stronger the mass loss rates, the faster the decline of the 
surface velocities. Let us add that polar winds may de- 
crease the quantity of angular momentum lost by stellar 
winds; however, as recalled in Sect. 2, we did not take ac- 
count of this effect due to the moderate initial velocities 
considered. 

How do these two processes, internal coupling and 
mass loss, depend on the initial metallicity ? For what 
concerns mass loss, it is now well established both theo- 
retically and observationally that the lower the metallicity, 
the lower the mass loss rates (Kudritzki & Puls l200ni Vink 
et al. '20015. This metallicity effect alone would produce a 
slower decline of the surface velocity at lower metallicity. 
Now what about the metallicity dependence of the inter- 
nal coupling mechanisms ? At lower metallicities, stars are 
more compact, and thus the densities in their outer layers 
at a given stage are higher. Since meridional circulation 
is more efficient than shear for the transport of angular 
momentum (Meynet & Maeder I2000|l , and since in the 
outer layers the driving effect for the meridional velocity 
is the "Gratton-Opik" term which is proportional to the 
inverse of the density, the outward transport of angular 
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Fig. 2. Evolution as a function of the actual mass of the 
rotation period, of the surface equatorial velocity and of 
the ratio of the angular velocity to the critical value during 
the WR stage of rotating stars. The long-dashed lines 
in the panels for the velocities show the evolution of the 
radius in solar units. Left: the WR phase of a star with 
an initial mass of 60 Mq with Vi^i — 300 km s^^ and Z — 
0.004. Right: for an initial mass of 60 Mq with Vini = 300 
km s-i and Z = 0.040. 



momentum in these layers will be less efficient at lower 
than at higher metallicity. This clearly reduces the inter- 
nal coupling, and weakens the value of the surface veloc- 
ity obtained at a given stage, all other things being equal. 
Thus when metallicity decreases, on one side, less angular 
momentum is removed by stellar winds, and on the other 
side, less angular momentum is brought up to the surface 
by the internal transport processes. Taken together these 
two effects tend to decrease the quantity of angular mo- 
mentum removed by stellar winds at low metallicity (see 
also Maeder & Meynet 12002 paper VII). 

For the evolution of the surface velocities, things are 
a little more complicated: at lower metallicity, the metal- 
licity dependence of mass loss rates favours higher values 
for the surface velocities at a given stage. On the other 
hand, the less efficient internal coupling favours lower val- 
ues for the surface velocities. The numerical models of the 
present work show that for the most massive stars (M > 30 
Mq), the effect of the metallicity dependence of the mass 
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loss rates dominates. This can be seen for instance in the 
case of the 40 M© tracks plotted in Fig. at high metal- 
licity [Z = 0.040) the surface velocity rapidly declines 
during the core H~burning phase, while at low metallicity 
(Z — 0.004) the star keeps a high surface velocity during 
most of this phase. 

The evolution of the ratio 17 /fic of surface angular ve- 
locity to the break-up or critical angular velocity is also 
very different at high and low metallicity. At a low metal- 
licity, the ratio fi/fic at the surface of the 60 M© model 
remains near the value of 0.5 during a great part of the 
Main Sequence phase, while at high metallicity the ratio 
continuously decreases. 

As was already emphasized in paper X, we again stress 
that any comparison between observed and predicted ro- 
tation for the large masses (M > 30 M©) is really much 
more a test bearing on the mass loss rates than a test 
of the internal coupling and evolution of rotation. For in- 
stance, the increase of Jl/ric at low metallicity was more 
pronounced in our previous low metallicity grids (papers 
VII and VIII) . This is because the mass loss rates adopted 
in these previous grids were lower by about a factor of two 
than the new values adopted here. This well illustrates 
both the uncertainties still pertaining to the mass loss 
rates and the sensitivity of the evolution of the surface ve- 
locity on this still not very well known physical ingredient 
of the stellar models. 

3.1. Evolution of the rotational velocities during the 
WR stages 

From Fig. ^ one can see that after the entry into the 
WR phase, the surface velocity can show rapid changes. 
These abrupt variations of the surface velocities are linked 
to rapid changes in the position of the star in the 
Hertzsprung-Russel diagram caused by surface compo- 
sition modifications induced by both rotational diffusion 
and mass loss by stellar winds. 

Fig. m shows specifically the evolution during the WR 
stages of the rotation period P (27r/f2), of the rotation 
velocities v at the equator and of the fraction ^ of the 
angular velocity to the critical angular velocity at the sur- 
face of a star model with an initial mass of 60 Mq and 
Vini — 300 km s~^ at low and high metallicities (left panels, 
Z = 0.004; right panels, Z = 0.040). One sees that at low 
metallicity, much greater and more rapid changes of the 
surface velocities are expected than at higher metallicity. 
Thus the encounter of the break-up limit during the WR 
phase may appear more often at lower metallicity (see the 
lower panels in Fig. |2J|. This comes from the fact that at 
low metallicity the whole WR phase occurs while the star 
is in the core He-burning phase, a phase during which our 
most massive star models experience strong changes in ra- 
dius. In contrast, at Z — 0.040, nearly the whole eWNL 
phase occurs during the core H-burning phase (for a def- 
inition of the eWNL phase see Sect. 5.1). One sees that 
the high mass loss rates undergone by the star makes the 



surface velocity to decrease regularly. Only when the star 
has lost its complete H-rich envelope (i.e. when the star 
enters the eWNE stage in terms of evolutionary models), 
does strong contraction of the surface layers, caused by 
the decrease of the opacity in the outer regions, produces 
a small increase of the surface velocity. We see that de- 
spite the strong increase of the velocity v and ^ do not 
change very much, since the radius is much smaller in the 
eWNE stage. All the transitions in the WR stage are fast, 
thus the transfer of angular momentum by meridional cir- 
culation is limited and the evolution of the rotation at the 
surface is dominated by the local conservation of angular 
momentum, which explains the fast and large decreases of 

p or n. 

One notes that in general the final surface velocity at 
low metallicity is also higher than at high metallicity, a 
fact which may be easily explained as a consequence of the 
lower mass loss rates at lower metallicity. From Tabled 
one sees also that at a given metallicity the velocities ob- 
tained at the end of the He-burning phase for the WR 
stars are in general smaller for the higher initial mass 
stars. This results from the fact that the higher the initial 
mass, the larger the amount of mass (and thus of angular 
momentum) removed by the stellar winds. 

The implications of the present models for the rotation 
rates of pulsars at birth and the possible progenitors of 
the collapsar models will be considered in a forthcoming 
paper. Let us just mention here that most of the WR 
stellar models have enough specific angular momentum 
in their core to be good candidates for collapsar models. 
On the other hand, if no angular momentum is lost during 
the core-collapse supernova explosion, these models would 
predict very fast rotational velocities for the neutron stars 
at birth. These results are in qualitative agreement with 
those obtained by Heger & Woosley H2004|l . and Woosley 
& Heger l|2003|l in the case of no magnetic field. 

4. Evolutionary tracks, lifetimes and final masses 

In Fig. |31 some evolutionary tracks of rotating models are 
shown. Note that during the WR phase the values of the 
effective temperature used to draw the tracks take account 
of the non-negligible optical thickness of the winds (see 
Sect. 1). At high metallicity and for the most massive 
stars, the entry into the WR phase occurs during the MS 
phase (see Sect. 5 below). Therefore the non-WR part of 
the track is quite short. For lower initial mass stars, rota- 
tion makes the tracks overluminous and more extended to- 
wards lower effective temperatures during the MS phase. 
This was also the case in previous works (see Heger & 
Laneer \MM Meynet & Maeder 001701 OIMIl . Note how- 
ever that for very fast rotation the star may be so effi- 
ciently mixed that it has a nearly homogeneous evolution 
and would present a very blue track during the MS phase 
(Maeder 1987). This is the case for the 500 km s^^ 60 M© 
stellar model at Z=0.004 in the present grid (see Table^. 
Let us recall also that at low metallicity, for those stars 
which do not already enter the WR phase during the MS 
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Fig. 3. Evolutionary tracks for 40 and 120 M© rotating models at different metallicities. The initial velocity is 300 km 
s~^. The light dotted lines correspond to the non-WR part of the tracks. The tracks during the WR phase are shown 
by heavy lines (continuous for the 120 Mq and dashed for the 40 Mq model). Symbols along the tracks are placed 
where the indicated surface hydrogen (Xsurf) and helium (Kjurf) abundances are reached. 



phase, the evolution towards the red supergiant stage is 
favoured when rotation is included in the stellar model 
(Maeder & Meynet EHIll ■ 

Table ^ presents some properties of the models. 
Columns 1 and 2 give the initial mass and the initial 
velocity v-mi respectively. The mean equatorial rotational 
velocity v during the MS phase is indicated in column 3. 
This quantity is defined as in paper V. The H-burning 
lifetimes tn, the lifetimes as an 0-type star on the MS 
to (we assumed that 0-type stars have an effective tem- 
perature higher than about 33 000 K), the masses M, the 
equatorial velocities v, the helium surface abundance Ys 
and the surface ratios (in mass fraction) N/C and N/0 at 
the end of the H-burning phase are given in columns 4 to 
10. Columns 11 to 17 present some characteristics of the 
stellar models at the end of the He-burning phase, tne is 
the He-burning lifetime. More details on the models are 
given in the electronic tables (see Sect. 2). 

From Table [l] one sees that for Z = 0.040 the MS 
lifetimes are increased by about 16-27% when the initial 



rotational velocity increases from Oto 300 kms-^ Similar 
increases were found at solar metallicity (see papers V and 
X). The He-burning lifetimes are increased by rotation 
for the 60 and 120 AIq models by respectively 12 and 
34%. This comes from the fact that the most massive stars 
already enter the WR phase during the core H-burning 
phase. The strong mass losses that they experienced at 
the end of the H-burning phase produce small He-cores 
at the beginning of the core He-burning phase. This tends 
to reduce the central temperatures and to increase the He- 
burning lifetimes. 

Rotation decreases the He-burning lifetime for the 25 
Mq by 8% at solar metallicity and by 17% at Z = 0.040. 
Contrary to what happens in the high mass star range, 
in the rotating 25 Mq stellar model the He-core at the 
end of the MS phase is significantly more massive than in 
the non-rotating model (by 42% a.t Z — 0.040) and thus 
the rotating track is much more luminous during most of 
the He-burning phase than the non-rotating one. This of 
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Table 1. Properties of the stellar models at the end of the H-burning phase and at the end of the He-burning phase. 
The masses are in solar mass, the velocities in km s^^, the lifetimes in million years and the abundances in mass 
fraction. 
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^ Model computed with anisotropic stellar winds and with an a— enhanced mixture of the heavy elements. 

^ Model computed with mass loss rates during the WR phase dependent on the mctallicity (Crowthcr et al. 2002). 



course tends to reduce the duration of the core He - burning 
phase. 

As indicated in Sect. 2, the present stellar models were 
computed with the mass loss rates of Vink et al. 1)2000(1 . 
This prescription shows a dramatic increase of the mass 
loss rates for OB Main-Sequence stars when they cross, 
from blue to red, the bistability limit which occurs at an 
effective temperature equal to ^25000 K. The increase 
is due to a drastic change in the ionization of the wind. 
To illustrate this effect, the evolution, during the Main- 
Sequence phase, of the mass of our stellar models is rep- 
resented in Fig. ^ as a function of Toft • When the star 
models cross the bistability limit, the tracks turn down 
very abruptly in the cases of the most massive stars, more 
smoothly in the cases of smaller initial mass stars. This 
results from the important increase of the mass loss rates, 
by about 1 dex, when the star crosses the bistability limit. 
As an example, the 25 Mq stellar model, which spends 
only 11% of its MS lifetime on the red side of the bistabil- 
ity limit, loses during this short phase more than 60% of 
the mass lost during the whole MS phase. For stars more 
massive than about 60 Mq and less massive than about 
10 Mq, there is no jump in the mass loss rates, since their 



tracks remain on one side of the bistability limit during 
the whole Main-Sequence phase. 

One can wonder to what extent the bistability effect 
represents an important feature favouring the entry of the 
stars in the WR regime. Let us just remark here that this 
effect alone is of little help in that respect. Indeed mod- 
els without rotation but taking account of this bistability 
effect give a poor fit of the observed variation of the num- 
ber ratio of WR to 0-type stars with the mctallicity (see 
Fig. ITUI) , while rotating models much better reproduce the 
observed data (see Sect. 6). 

From Tabled comparing the final masses obtained at 
Z = 0.040 for the rotating and the non-rotating model, 
one can see that rotation produces in general smaller final 
masses. Fig.|Slshows the final masses obtained for different 
initial masses and metallicities. K\ Z — 0.040, most of 
the stars in the mass range considered here end their life 
with masses below 10 Mq. All stars with masses above 
about 50 Mq reach a final mass between 5 and 7.5 Mq. In 
contrast, at a mctallicity one order of magnitude smaller, 
the final masses of stars are significantly higher, being in 
the range between 17 and 29 Mq for stars with initial 
masses above 60 Mq . This well illustrates the effect of the 
mctallicity dependence of the mass loss rates, an effect 



G. Meynet, A. Maeder: Stellar evolution with rotation XI: 



9 




Fig. 4. Evolution of the mass of the stars as a function of 
the effective temperature for different initial mass models 
(from 9 to 40 Mq) at various metallicities during the Main- 
Sequence phase (uini = 300 km s~^). The initial mass of 
the stars is the ordinate of the hottest point of each track 
(left point). Models at Z=0.004, 0.008, 0.020 and 0.040 are 
shown with dotted, short-dashed, continuous and long- 
dashed lines respectively. One can note the effect of the 
crossing of the bistability limit around 25000 K (see text). 

which is qualitatively similar to that found in non-rotating 
stellar models fMaeder [l991|l . 

5. Rotating models of Wolf-Rayet stars 

5.1. Effects of rotation and metallicity on the 

evolutionary scenarios leading to the formation of 
WR stars 

The effects of rotation on the evolution of massive sin- 
gle stars into the Wolf-Rayet phase have been discussed 
by Maeder H1987|l . Fliegner and Langer l|1995|l . Maeder 
& Mevnet 12000^ and Meynet ipUUjl . These studies were 
mainly based on solar metallicity models. Here we extend 
these discussions to lower and higher than solar metallic- 
ities. 

We had to choose a set of criteria to decide when a 
stellar model enters the WR phase. Ideally, of course, the 
physics of the models should determine when the star is 
a WR star. However our poor knowledge of the physics 
involved, as well as the complexity of models coupling the 
stellar interiors to the winds, are such that this approach 
is not yet possible. Instead, it is necessary to adopt some 
empirical criteria for deciding when a star becomes a WR 
star. In this work the star enters the WR phase when 
two conditions are fullfilled: 1) log Toff > 4.0 (note that 



Fig. 5. Relations between the final and the initial mass 
for rotating stellar models at various metallicities. The 
line with slope one, labeled M = 0, corresponds to the 
case without mass loss. 

after the entry into the WR phase the star may excep- 
tionally spend very short periods in cooler regions of the 
HR diagram. In general however, the whole WR phase is 
spent at effective temperatures well above log Toff — 4.0); 
2) the mass fraction of hydrogen at the surface Xs is in- 
ferior to 0.4. Reasonable changes to these values (for in- 
stance adopting < 0.3 instead of 0.4) do not affect 
the results significantly. Let us recall that the effective 
temperature considered here is an effective temperature 
appropriately averaged over the surface (see Meynet & 
Maeder I1997II . We consider the WR star to be of the 
eWNL type when the mass fraction of hydrogen at the 
surface is superior to 10~^, (we adopt here the denom- 
ination proposed by Foellmi et al. I2003bl for the WNL 
and WNE phase, based on evolutionary criteria instead of 
spectroscopic ones, hence the small "e" before the name). 
The eWNE phase begins at the end of the eWNL phase, 
while the transition WN/WC phase is considered to begin 
when the mass fraction of carbon at the surface becomes 
superior to 10% of the mass fraction of nitrogen. Both the 
eWNE and the transition WN/WC phases end when the 
WC/WO phase begins, which we assume to begin when 
the mass fraction of nitrogen at the surface becomes less 
than 10% of the mass fraction of carbon. We shall not dis- 
tinguish here between WC and WO stars. Also, changes 
of the numerical values of the limits have very little con- 
sequence since the transitions are fast. 

Internal mixing favours the entry into the WR phase in 
two ways, firstly by allowing chemical species produced in 
the core to diffuse in the radiative envelope and, secondly, 
by making the mass of the convective core larger (see also 
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Fig. 6. Evolution as a function of time of the total mass 
Mxot and of the mass of the convective cores Mcc during 
the H- and He-burning phases, for a 60 Mq stellar model 
with and without rotation a,t Z = 0.040. The type of the 
WR star at a given age is given in the upper part of the 
figure. 

Meynet & Maeder FiUOUfl . In the non-rotating model, mass 
loss by stellar winds is the key physical ingredient which 
allows internal chemical products to appear at the surface 
and thus form a WR star. The star becomes a WR star 
only when sufficiently deep layers are uncovered. 

In Fig. ini the evolution of the structures for 60 Mq 
models with and without rotation at twice the solar metal- 
licity are shown. In the case of the 60 Mq model, the most 
striking differences between the non-rotating and the ro- 
tating model are the following: 

— Inclusion of rotation in the models allows an earlier 
entry into the WR phase. Typically here in the non- 
rotating model the star becomes a WR star when the 
actual mass is 34.3 Mq and the mass fraction of H at 
the centre X,, is 0.04, while in the rotating model it 
enters the WR phase when the mass is equal to 44.9 
Mq and Xc is equal to 0.24. 

— As a consequence the WR lifetime as well as the dura- 
tion of the eWNL phase will be increased by rotation. 
Also the actual luminosity at the entry of the WR 
phase will be higher for the rotating models. 



— The duration of the WC phase is reduced. 

Qualitatively similar effects of rotation were obtained at 
solar metallicity (cf. paper X) except for the effect of rota- 
tion on the WC phase: at solar metallicity and for the 60 
M0 stellar model the WC phase is slightly enhanced by 
rotation (see Table 2 in paper X). At Z = 0.004, the WC 
phase becomes longer when the initial rotation increases 
(see Table |3J . In addition to the above effects, rotation 
decreases the minimum initial mass of single stars going 
through a WR phase (see Fig. |HJ| and the duration of the 
transition WN/WC phase is much longer in the rotating 
models than in the non-rotating ones, (see Table [SJ. 

As said in Sect. 2, the importance of these effects de- 
pends on the values adopted for the initial velocity and 
on various physical ingredients of the stellar models as 
for instance the mass loss rate and/or the overshooting. 
Some illustration of the effects of different values of rota- 
tion for 60 M0 models at solar metallicity can be found in 
Meynet ffOW) and in paper X. Results for a fast rotating 
60 Mq stellar model at Z=0.004 (vini = 500 km s'^) are 
presented in Tables ^ and 13 In that last example, we see 
that increasing the initial velocity from 300 to 500 km s~^ 
more than doubles the WR lifetime. This enhancement 
results from the more efficient internal mixing at high ro- 
tation and not to the increase of the mass loss rate due to 
rotation. 

From the present rotating models one can derive two 
interesting limiting masses. The first, MqwRj is the mini- 
mum initial mass of a single star entering the WR phase 
during the MS phase. The second, Mwr, is the minimum 
initial mass of a single star entering the WR phase at 
any point in the course of its lifetime. These two limit- 
ing masses define the mass ranges of three evolutionary 
scenarios for the massive stars: 

— For M > MowR, the stars will avoid the Luminous 
Blue Variable stage after the MS phase. In this case, 
they will go through the following phases: 0-eWNL- 
eWNE-WC/WO. Note that here we assume that once 
the star has entered the WR regime, it remains a WR 
star for the rest of its lifetime. Some stars however 
may evolve in cooler regions of the HR diagram af- 
ter they have entered the WR phase and might thus 
encounter the Humphreys-Davidson limit. These star 
models would present characteristics similar to LBV 
stars and could thus belong to this category. 

— For Mwr < M < Mqwr, after the MS phase, the 
star will evolve into the cooler part of the HR dia- 
gram, where it may encounter the OF-limit (Maeder 
& Meynet ;200(Q or become a Red Supergiant. In that 
case, one would have 0-LBV or RSG-eWNL-eWNE- 
WC/WO. Evolution may not necessarily proceed up 
to the WC/WO stage, it may stop at the eWNE or 
eWNL stage. 

— For M < Mwr, after the core H-burning phase the O- 
type star will become a supergiant, but it will never 
enter the WR phase. The blue or red nature of the su- 
pergiant depends among other parameters on rotation 
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(Maeder & Mevnet l!jUl)l|l . The evolutionary sequence 
in that case will be 0-RSG/BSG. Blue loops may be 
present especially in the lower mass star range. 

In order to estimate Mqwr and Mwr we interpolated 
between corresponding evolutionary stages of the tracks 
using the logarithm of the initial mass as the interpolat- 
ing factor. We obtain the mass limits indicated in Table[21 
We can note that, at a given metallicity, the mass limits 
are lower when the rotation rate is higher, as expected 
from the considerations mentioned above. As we shall see 
in Sect. 6, these results on the mass limits have interest- 
ing consequences for the expected number of WR stars at 
different metallicity, as well as for the predicted number of 
type Ib/Ic supernovae. One notes also from Table El that 
the mass limits are lower at higher metallicities, as would 
be the case for models without rotation. This is due to 
the greater mass loss rates experienced by stars at higher 
metallicity. 

Table 2. Mass limits for the different evolutionary sce- 
narios (see text). 

Z Uini Mwr Mqwr 

[km s-i] [A/q] [Mq] 



0.040 





29 


42 


0.040 


300 


21 


39 


0.020 





37 


62 


0.020 


300 


22 


45 


0.008 


300 


25 


69 


0.004 


300 


32 


75 



The mass range for which the WR phase is preceded 
by a LBV and/or a RSG phase shifts to higher values at 
lower metallicity, it also extends over a larger mass inter- 
val. Typically, at Z = 0.040 rotating models predict that 
such a scenario for single massive stars occurs in the mass 
range between 21 and 39 Mq, while it occurs in the mass 
range between 32 and 75 Mq at Z = 0.004. Since qual- 
itatively a same trend is predicted by the non-rotating 
models, this result would imply that the upper luminos- 
ity of LBV stars should be lower at higher metallicity. 
However one should be cautious here. First, one cannot 
discard the possibility that a star enters the WR phase and 
later undergoes a shell ejection similar to a LBV outburst 
due for instance to reaching the riF-limit (see Maeder & 
Mevnet r2000a|l . Moreover the observed stellar populations 
are a mixture of stars of different initial velocities, whose 
initial distributions may depend on the metallicity. This 
would contribute to the blurring of the schematic picture 
just described above. 

As already stated in Sect. 2, the results shown is 
Table El are sensitive to various physical ingredients. For 
instance, an enhancement of the core sizes and/or of the 



initial velocity and/or of the mass loss rates would favour 
the WR formation and thus the lowering of these mass lim- 
its. Thus the numbers shown in TableOmay undergo some 
changes in the future when improvements in these physical 
ingredients will be reached. However, at least qualitatively, 
the trend they show, namely the fact that rotation lowers 
the mass limits, will likely persist. This is well illustrated 
by the fact that a similar trend was already present in 
our paper V models which were computed with different 
prescriptions for the mass loss rates, the overshooting and 
the rotational mixing (see Sect. 2). 

For a given initial mass and an initial rotational veloc- 
ity, the greater the metallicity the larger the luminosity 
range spanned during the WR phase, and a lower mini- 
mum luminosity is reached. This can be seen from Fig.O 
Interestingly we can see that at a given metallicity, the 
higher initial mass stars are in general the progenitors of 
the less luminous WC stars. This is mainly due to the fact 
that the mass loss rates increase with the luminosity, en- 
abling the most massive stars to enter the WR phase at 
an earlier stage and thus to lose mass at high rates dur- 
ing a larger portion of their lifetime. Thus these stars end 
with small final masses and hence with small luminosities 
since at this stage they are following a well known mass 
hmiinosity relation (Schaerer & Maeder [1992(1 . 

5.2. The Wolf-Rayet lifetimes 

The WR lifetimes of the present rotating models for the 
four metallicities are plotted as a function of the initial 
mass in Fig. |3 As was the case with the non-rotating 
models fMaeder lTMTl Maeder & Mevnet [TUM|l . the metal- 
licity dependence of the mass loss rates is responsible for 
two features: 1) for a given initial mass and velocity the 
WR lifetimes are greater at higher metallicities. Typically 
at Z=0.040 and for M > 60 Mq the WR hfetime is of the 
order of 2 Myr, while at the metallicity of the SMC the 
WR lifetimes in this mass range are between 0.4-0.8 Myr; 
2) the minimum mass for a single star to evolve into the 
WR phase is lower at higher metallicity. 

It is interesting to compare the WR lifetimes obtained 
dX Z — 0.040 and 0.004 from different sets of models (see 
Fig.jSl). First, one can notice that the present non-rotating 
models with the updated mass loss rates give very similar 
results to those obtained with the "normal" mass loss rates 
by Meynet et al. I|1994l see how the black squares give a 
good fit to the Z — 0.040 dotted curve). As was discussed 
in paper X, rotation has a similar effect as an enhance- 
ment of the mass loss rates on the WR lifetimes. Namely, 
for a given initial mass it increases the WR lifetime and 
also lowers the minimum initial mass of single stars going 
through a WR phase. Note however that quantitatively 
the results obtained from rotating models may be differ- 
ent from the one obtained from the enhanced mass loss 
rate models. At low metallicity, the curve obtained from 
the rotating models is always above the one obtained from 
the enhanced mass loss rate models. More importantly, 
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Fig. 7. Lifetimes of Wolf-Rayet stars from various initial 
masses at four different metallicities. All the models begin 
their evolution with Umi — 300 km s~^ on the ZAMS, the 
corresponding time-averaged velocities during the MS O- 
type star phase are given in Tabled 
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Fig. 8. Lifetimes of Wolf-Rayet stars from various initial 
masses for Z = 0.040 and 0.004 from different sets of mod- 
els: continuous lines are for the present rotating models 
with v\ai = 300 km s^^, the dotted and dashed lines are for 
the non-rotating stellar models with "normal" and "en- 
hanced" mass loss rates computed by Meynet et al. p994|l . 
The black squares indicate the lifetimes obtained for the 
non-rotating 25, 60 and 120 Mq models obtained in the 
present work at Z = 0.040. 



rotation seems to be much more efficient in lowering the 
value of MwR. This is due to the fact that rotation helps 
in forming WR stars through internal mixing, a process 
which, in contrast to mass loss rates, is important for all 
massive stars, not only the most massive ones. At high 
metallicity, the situation is somewhat different, because 
when metallicity increases the mass loss rate becomes the 
dominant parameter. 



In Fig. El the duration of the different WR subphases 
is plotted as a function of the initial mass for the vari- 
ous metallicities considered in this work. Only the results 
with rotation are plotted. The durations of these differ- 
ent subphases are given in Table 01 The greatest part of 
the WR lifetime is spent in the eWNL phase. Note that, 
as explained above, rotation increases the duration of this 
phase. For identical initial rotational velocities, the dura- 
tion of the eWNL phase is greater at higher metallicity. 
This results from the metallicity dependence of the mass 
loss rates. At higher metallicity, the higher mass loss rates 
by stellar winds enable the star to enter the WR phase at 
an earlier stage. The eWNE phase is also longer at higher 
metallicity (as is also the case for non-rotating models). 
The WC phase keeps more or less the same duration for 
all the metallicities in the higher mass star range. In the 
lower mass star range the WC phase is longer at higher 
metallicity as a result of the shift towards a lower value of 
the minimum initial mass of single stars needed to become 
a Wolf-Rayet star. 

As already explained in paper X, in the rotating stel- 
lar models a new phase of modest, but non-negligible du- 
ration, appears: the so-called transition WN/WC phase 
(see the values of ^wn/wc given in Table This phase 
is characterized by the simultaneous presence at the sur- 
face of both H- and He-burning products. The reason for 
this is the shallower chemical gradients which build up in- 
side the rotating models. These shallower gradients inside 
the stars also produce a smoother evolution of the surface 
abundances as a function of time (see Fig. I14|l . As ex- 
plained in paper X, for a transition WN/WC phase to oc- 
cur, it is necessary to have — for a sufficiently long period — 
both an He-burning core and a CNO-enriched envelope. 
In general, in the highest mass stars, mass loss removes the 
CNO-enriched envelope too rapidly to allow a long tran- 
sition WN/WC phase to occur. In the low mass range, 
the time spent in the WR phase is too short and the H- 
rich envelope too extended to allow He-burning products 
to diffuse up to the surface. Consequently, a significant 
transition WN/WC phase only appears in an intermedi- 
ate mass range between ~ 30 and 60 M© for Wini = 300 
km s~^. In general the duration of this phase also increases 
significantly when mass loss rates depending on the metal- 
licity are used during the WR phase (see Sect 7). 

6. Comparison with the observed properties of 
Wolf-Rayet star populations 

Following Maeder H1991(l and Maeder & Meynet 1] 1994(1 . 
one can easily estimate the theoretical number ratio of WR 
to 0-type stars in a region of constant star formation. 
This ratio is simply given by the ratio of the averaged 
lifetimes of a WR star to that of an OV-type star. The 
averaged lifetimes are weighted means of the lifetimes over 
the initial mass function (IMF). Assuming a Salpeter IMF 
slope (diV/dM oc M~'^^+^\ with x = 1.35), considering 
the 0-type and WR star lifetimes given in Tables d and 
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Fig. 9. Variation of the durations of the WR subphases as a function of the initial mass at various metalhcities. All 
the models begin their evolution with Vi^i = 300 km s^^. 



Inland the mass limits A/wr given in Table we obtain 
the predicted ratios given at the bottom of Table |2| 

Let us emphasize here that in the present work we as- 
sume that the Vmi = 300 km s^^ stellar models are well 
representative of the behaviour of the majority of the OB 
stars. In future studies, when information on the distri- 
bution of the initial rotation velocities will be available 
at different metalhcities, it will certainly be interesting to 
convolve the theoretical results obtained for various initial 
velocities with these observed velocity distributions. 

Looking at the results presented in Table |31 one can 
note that rotating models predict that 

— The fraction of WR stars with respect to 0-type stars 
increases with the metallicity. 

— The fraction of eWNE with respect to WR stars in- 
creases with the metallicity, as does the fraction of WC 

At Z — 0.004, the value for the minimum initial mass of 
O-type stars is taken from the rotating models of Maeder & 
Meynet ll^IMTl . At Z = 0.008 we use the value obtained from 
the "enhanced" mass loss rate models of Meynet et al. 119941 
see Table 2 in Maeder & Mevnet 119911 . 



stars (although the ratio seems to saturate around 0.25 
for Z > 0.020). 

— The fraction of eWNL decreases with the metallicity 
again reaching a saturation level around 0.65 for Z > 
0.020. 

— The fraction of stars in the transition WN/WC phase 
has a non-monotonic behaviour, increasing from 0.01 
at Z = 0.004 to 0.07 at Z = 0.008, and decreasing to 
values between 0.04-0.05 for Z > 0.020. 

Comparisons with observed number ratios are shown 
in Figs. 1101 and llll The values for the non-rotating models 
(both those computed with the present physics and those 
computed by Maeder & Mevnet fl994(l are well below the 
observed values. The ratios predicted by the models with 
rotation are in much better agreement with the observa- 
tions. In particular the ratios predicted by the rotating 
models are in quite good agreement with the observed 
ones at low metallicity. They give predictions that are well 
in line with the results obtained by Foellmi et al. H2003al 
I2003b|l . Let us recall that these authors concluded that, 
in contrast to previous expectations, the binary channel 
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Table 3. Lifetimes in Myr of WR stars of different ini- 
tial masses and metallicities. The durations of various WR 
subphases are also indicated as well as theoretical predic- 
tions for various number ratios. 



M 


^ini 




^ C VV In Li 


r--i'\A7'T\.T "FT 
^ C VV IN 


-■ WN / WC 










Z= 


=0.004 






120 


300 


0.8298 


0.5105 


0.0066 


0.0089 


0.3127 


60 


300 


0.3704 


0.3126 


0.0039 


0.0113 


0.0539 


60 


500 


0.8617 


0.4617 


0.1422 


0.1089 


0.2578 


40 


300 


0.3057 


0.3057 


0.0000 


0.0000 


0.0000 








/. 


=0.008 






120 


300 


1.1107 


0.7354 


0.0282 


0.0167 


0.3471 


60 


300 


0.4315 


0.2113 


0.0088 


0.0186 


0.2114 


40 


300 


0.4677 


0.4168 


0.0238 


0.0372 


0.0271 


30 


300 


0.3421 


0.2971 


0.0299 


0.0528 


0.0151 








Z= 


=0.040 






120 





1.1829 


0.7646 


0.1276 


0.0034 


0.2907 


120 


300 


1.8956 


1.4398 


0.2004 


0.0048 


0.2554 


120^ 


300 


2.3190 


1.7200 


0.4339 


0.0382 


0.1651 


85 


300 


1.8955 


1.4626 


0.1641 


0.0980 


0.2688 


85 1 


300 


2.2710 


1.6190 


0.3087 


0.1723 


0.3433 


60 





0.5592 


0.2083 


0.0552 


0.0009 


0.2957 


60 


300 


1.8083 


1.3294 


0.2640 


0.0083 


0.2149 


60^ 


300 


2.1382 


1.5450 


0.4515 


0.0940 


0.1417 


40 


300 


0.5960 


0.2857 


0.0537 


0.0291 


0.2566 


40^ 


300 


0.5710 


0.2221 


0.0433 


0.0231 


0.3056 


25 


300 


0.4493 


0.2091 


0.0706 


0.0602 


0.1696 








Predicted number ratios 
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cWNL 


cWNE 


WN/WC 


WC 


WC 2 




WR 


WR 


WR 


WR 


WN 








z= 


=0.004 






300 


0.02 


0.83 


0.00 


0.01 


0.15 


0.18 








z= 


=0.008 






300 


0.05 


0.74 


0.05 


0.07 


0.22 


0.28 








z= 


0.020^ 









0.02 


0.35 


0.16 


0.00 


0.49 


0.97 


300 


0.07 


0.66 


0.05 


0.04 


0.25 


0.35 








z= 


=0.040 









0.06 


0.46 


0.10 


0.00 


0.44 


0.78 


300 


0.16 


0.63 


0.12 


0.05 


0.24 


0.32 


300^ 


0.18 


0.59 


0.15 


0.06 


0.26 


0.36 



^ Model eomputed with mass loss rates d uring the WR pliase dependent 
on the mctalUcity (Crowther et al. 12001^1 . 

2 WC 



cWNL+cWNE ■ 

Values obtained in paper X. 



for WR star formation does not seem to be a favoured 
scenario at low metallicity. Indeed, as already mentioned 
in the introduction, the fraction of binaries among the 
WR stars of the Small (40%) and Large (30%) Magellanic 
Cloud appears to be similar to that found in our Galaxy. 
Moreover, not all WR stars in a binary system must nec- 
essarily owe their existence to a Roche Lobe Overflow, 
since the two stars in a binary system may be sufhciently 
distant from each other to prevent that such a process 
occurs. Thus models from single massive stars should be 
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Fig. 10. Variation of the number ratios of Wolf-Rayet 
stars to 0-type stars as a function of the metallicity. The 
observed points are taken from Maeder & Meynet H1994|l . 
The dotted line shows the predictions of the models of 
Meynet et al. 1)1994(1 with normal mass loss rates. The 
continuous and the dashed lines show the predictions of 
the present rotating and non-rotating stellar models re- 
spectively. The black pentagon shows the ratio predicted 
by Z=0.040 models computed with the metallicity depen- 
dence of the mass loss rates during the WR phase proposed 
by Crowther et al. ijMT^ . 
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Fig. 11. Variation of the number ratios of WN to WC 
stars as a function of metallicity. The black circles are ob- 
served points taken from Massey & Johnson 1)20011 and see 
references therein), except for the SMC (Massey & Duffy 
IWr|> . for NGC 300 (Schild et al. IMT^ and for ICIO, 



for which we show the estimate from Massey & Holmes 
(2002 ) . The continuous and dotted lines show the pre- 
dictions of the present rotating and non-rotating stellar 
models respectively. The black pentagon shows the ratio 
predicted by Z=0.040 models computed with the metal- 
licity dependence of the mass loss rates during the WR 
phase proposed by Crowther et al. 1)2002)1 . 
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Fig. 12. Variation of the number ratios of type Ib/Ic su- 
pernovae to type II supernovae. The crosses with the error 
bars correspond to the values deduced from observations 
by Prantzos & Boissier Ij2003|l . The dotted line is an ana- 
lytical fit proposed by these authors. The continuous and 
dashed line show the predictions of the present rotating 
and non-rotating stellar models. 



able to reproduce more than half the number of observed 
WR stars in the Magellanic Clouds, a condition quite well 
fulfilled by the rotating models and not realized by the 
non-rotating ones. At high metallicity, the agreement is 
less good, although still much better than the one obtained 
from non-rotating models. Let us note that the samples 
at high metallicity do not have the same level of com- 
pleteness as for the Magellanic Clouds. In the Clouds, the 
distance and metallicity are relatively well known quan- 
tities, moreover the internal extinction is relatively weak 
and thus does not constitute a severe barrier. The central 
regions of galaxies and of the Milky Way in particular do 
not offer such good conditions, and completeness problems 
are likely much more important than in the Magellanic 
Clouds. 

In Fig.ini the predicted values for the WC/WN num- 
ber ratio are compared with observations in regions sup- 
posed to have undergone a constant star formation rate. 
The observations show that the WC/WN ratio increases 
with the metallicity along a relatively well defined relation. 
The observed point for the solar neighborhood is however 
well above the general trend. According to Massey 1)200311 
this may result from an underestimate of the number of 
WN stars. 

At low metallicity, models with rotation are well within 
the general observed trend. At solar metallicity they are 
just at the inferior limit and at twice the solar metallic- 
ity they are below the extrapolated observed trend. The 
non-rotating models predict (at least for Z > 0.020) sig- 
nificantly higher WC/WN ratios. If the value predicted for 
the solar metallicity lies well above the real observed value 
at this metallicity (see Massev 120031) . the one at twice the 



solar metallicity appears to be in good agreement with the 
extrapolated observed trend. 

Taken at face value, it seems that rotating models 
overestimate the number of WN stars at high metallic- 
ity. Apart from invoking completeness problems in the 
observed sample, how might the comparisons between the 
ratios predicted by the rotating models and the observed 
ones be improved ? We see two possible solutions: 1) From 
Fig. 111! it clearly appears that non-rotating stars, or more 
reasonably, stars rotating with relatively small velocities 
at high metallicity, would make the predicted points re- 
main in the general observed trend. Slow rotators undergo 
little rotational mixing, thus they keep steep gradients of 
chemical composition in their interior. The fraction of the 
mass of the star presenting abundances characteristic of 
those of WN stars is therefore smaller than in a rotating 
star. This tends to decrease the WN lifetime. 2) The use 
of higher mass loss rates during the WN phase will tend to 
remove the H-rich envelope more rapidly, decreasing the 
WN lifetime and favouring an early entry into the WC 
phase. 

The first solution is not reasonable for at least two 
reasons. Firstly, Fig. ^| shows that slowly rotating mod- 
els are unable to account for the high fraction of WR stars 
observed at high metallicity. Secondly the observed num- 
ber ratio of type Ib/Ic supernovae with respect to type 
II supernovae would also be severely underestimated by 
non-rotating models (Prantzos & Boissier 120031 see also 
Sect. 6.1 below). The second solution, enhanced mass loss 
rates during the WN phase, appears more promising, al- 
though the effects of such an enhanced mass loss rate may 
be different according to the stage at which the star enters 
the WR phase. The most massive stars already enter the 
WN phase during the MS phase. An enhancement of the 
mass loss rates in that case can have the opposite effect of 
the one desired, namely it may lengthen the WN phase. 
Indeed a higher mass loss produces a more important re- 
duction of the convective H-burning core, accompanied by 
a severe reduction of the central temperature. This tends 
to make the end of the H-burning phase longer, and thus 
to eventually increase the duration of the WN phase. Only 
an enhancement of the mass loss rates for those WN stars 
beyond the H-burning phase would reduce the WN phase. 

Recently Crowther et al. H2002|l proposed a metallicity 
dependence for the mass loss rates of WR stars compara- 
ble to that predicted for 0-type stars, i.e. with (Z/Z©)^/^. 
We have recomputed the Z = 0.040 models between 40 
and 120 Mq with this new prescription. The lifetimes of 
the different WR subphases are given in Table 13 We note 
that for the 120 Mq star, which enters the WR phase 
during the core H-burning phase, the stronger mass loss 
rates given by the new prescription lead to an increase 
of the WN phase and a decrease of the WC phase. The 
reason is just the one invoked above. In the case of the 
40 Mq model, where only 17% of the WN phase is spent 
in the core H-burning phase, the metallicity dependent 
mass loss rates during the WR phase give a shorter WN 
phase and a slightly longer WC phase. It is likely that in 
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the case of the 25 Mq model, which spends its whole WR 
phase in the post core H-burning phase, the reduction 
of the WN phase and the enlargement of the WC phase 
would be more important. Supposing that for this 25 Mq 
stellar model the total WR lifetime is not much changed 
by the use of metallicity dependent mass loss rates dur- 
ing the WR phase, one can estimate the durations of the 
WN and WC phases one would obtain by simply assum- 
ing that the new iwN is equal to that given in Table 01 
divided by (0.04/0.02)^/^ = 1.4 and the new tYiC is equal 
to iwR — ^WN- If one docs this one finds that the new twN 
and twc are 0.20 and 0.25 Myr respectively. With these 
values, one obtains the predicted number ratios given in 
the last line of TableOl One sees that the use of metallicity 
dependent mass loss rates during the WR phase does not 
much increase the predicted value for the WC/WN ratio 
at high metallicity (see also Fig. . 

Thus in view of the above discussion, the two possibil- 
ities which presently appear as the most probable for rec- 
onciling the predicted WC/WN ratio with the observed 
ones are 1) that the observed WC/WN ratio is overes- 
timated as seems to be the case at least for the Milky 
Way fMassev l2003|) : 2) that the mass loss rates during 
the WN phase after the core H-burning phase are higher 
than presently assumed. One may also wonder what is the 
role played by binaries in this context. 

6.1. The type Ib/lc supernovae 

Current wisdom associates the supernovae of type Ib/Ic 
with the explosion of WR stars, the H-rich envelope of 
which has been completely removed either by stellar winds 
and/or by mass transfer through Roche Lobe overflow in 
a close binary system. If we concentrate on the case of 
single star models, theory predicts that the fraction of 
type Ib/Ic supernovae with respect to type II supernovae 
should be higher at higher metallicity. The reason is the 
same as the one invoked to explain the increasing num- 
ber ratio of WR to 0-type stars with the metallicity Z, 
namely the growth of the mass loss rates with Z. Until 
recently very little observational evidence has been found 
confirming this predicted behaviour. The situation began 
to change with the work by Prantzos & Boissier (2003). 
These authors have derived from published data the ob- 
served number ratios of type Ib/Ic supernovae to type 
II supernovae for different metallicities. The regions con- 
sidered are regions of constant star formation rate. Their 
results are plotted in Fig. E| Using the predicted ratios at 
solar metallicity deduced from our rotating models (paper 
X), they concluded that rotating models were much bet- 
ter able than non-rotating ones to reproduce the observed 
ratio at solar metallicity. 

We complement here the discussion of Prantzos & 
Boissier (12003|l by comparing theory to the observation 
at other metallicities. Theoretical ratios are simply given 
by the ratio A^ib/ic/-^//i where iVib/ic is equal to the inte- 
gration of the IMF over the range of Mwne to 120 Mq, 
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Fig. 13. Evolution during the MS phase of the N/C ratios 
(in number) at the surface of rotating stellar models as a 
function of the effective temperature. The differences in 
N/C ratios are given with respect to the initial values. 

where Mvvne is the minimum initial mass of the stars that 
end their lives as eWNE or WC stars, and where Nu is 
the integration of the IMF between 8 and Mwne- To de- 
termine the values of Mwne we used the same technique 
as the one used to determine Mwr and Mqwr- We ob- 
tained 20.6, 22, 25 and 52 for Z ^ 0.040, 0.020, 0.008, 
0.004 respectively. Looking at Fig. ^1 it clearly appears 
that the conclusion of Prantzos & Boissier H2003|l is con- 
firmed by the comparisons at low and high metallicities. 
Rotating models give a much better fit to the observed 
data than non-rotating models. This comparison can be 
viewed as a check of the lower initial mass limit Mwne 
of the stars evolving into a WR phase without hydrogen, 
while the comparison between the observed and predicted 
number ratio of WR to 0-type stars involves not only the 
value of the minimum initial mass of stars evolving into 
the WR phase but also the durations of the WR phase. 
In that respect the comparison with the supernovae ratios 
is a more direct check of the correctness of the value of 
Mwne which at high metallicity appears rather close to 
Mwr. 

7. Evolution of the surface abundances 

The effects of rotation on the surface abundance of massive 
stars at solar composition have been discussed by Heger 
& Langer ipTOjl . Meynet & Maeder 1^017011^017511 . One of 
the main features is the enhancement of the N/C ratios at 
the surface during the MS phase. In papers VII and VIII 
we obtained that, for a given initial mass and velocity, the 
surface enrichment induced by rotation is higher at lower 
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Fig. 14. Evolution as a function of the actual mass of the abundances (in mass fraction) of different elements at the 
surface of rotating 40 Mq stellar models at various metallicities. 
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metallicity. This is a consequence of the fact that at lower 
metallicity stars are more compact and have steeper inter- 
nal gradients of 17 (see paper VII and VIII). This implies 
strong shear mixing of the chemical elements and explains 
why in Fig. E| the N/C ratio at the surface of the 9 Mq 
model at Z = 0.004 is higher than for the corresponding 
model at Z = 0.020. 

For the more massive stars the situation is different. 
Indeed one can see in Fig.^Jthat the N/C ratio at the end 
of the MS phase at the surface of the rotating 40 Mq at 
Z=0.040 is higher than the N/C ratio obtained at the same 
stage at the surface of the 40 Mq at Z=0.004. Why is the 
behaviour of the most massive stars so different from that 
of the less massive ones ? The reason is that, in the high 
mass star range, stellar winds become the dominant effect. 
High mass losses facilitate the surface enrichment by un- 
covering deep layers whose chemical composition has been 
changed by nuclear processing and/or rotational mixing. 

Fig. 1141 shows the evolution of the surface abundances 
in rotating models with an initial mass of 40 M©. Before 
commenting on the effects of the metallicity, let us briefly 
recall the effects of rotation. In rotating models one can 
note the following differences with respect to the non- 
rotating ones: 

— In rotating models the progressive changes of the abun- 
dances of CNO elements from the initial cosmic values 
to the values of the nuclear equilibrium of the CNO 
cycle are smoother. This is due to rotational mixing 
which smoothes internal chemical gradients. 

— Due to rotational mixing, the change of abundances 
also occurs much earlier in the peeling-off process. 
This is also true for the changes of H and He. We 
note that the nuclear equilibrium CNO values are es- 
sentially model independent as already stressed a long 
time ago (Smith & Maeder ri991|l . This is true whether 
H is still present or not. 

— In the rotating case, the transitions between the eWNE 
phase and the WC stage are smoother, so that there 
are some stars observed in the transition state which 
correspond to the so-called WN/WC stars, (Conti & 
MasseyHnnni Crowther et al. ll995;i . These stars show 
simultaneously some ^^C and '^''N. They could also 
have some ^^Ne excess. Since the attribution of spec- 
tral types is a complex matter, it may even be that 
some of the stars in the transition stage are given a 
spectral type WNE or WC, thus we might well have a 
situation where a WN star would have some ^^Ne ex- 
cess or a WC star would still have some ^"^N present. 

— At the entry into the WC phase, the ^^C and ^^O abun- 
dances are lower in rotating models, and the abun- 
dance of He is higher. 

— Also, the fraction of the WC phase spent with lower 
C/He and 0/Hc ratios is longer in models with rota- 
tion. 

— However, rotation does not affect the high level of the 
^^Ne abundance during the WC phase. This is a conse- 
quence of the fact that most of this ^^Ne results from 



the transformation of the ^''N produced by the CNO 
cycle in the previous H-burning core. The value of the 
^''N and therefore that of the ^^Ne is fixed by the char- 
acteristics of the CNO at equilibrium, which in turn 
depends on the nuclear physics and not on the pecu- 
larities of the stellar models. It is interesting to men- 
tion here that the high overabundance of ^^Ne at the 
surface of the WC star predicted by the models is well 
confirmed by the observations f Willis [T999I Dessart et 

ai.Ennni). 

Looking at Fig. one can see how the evolution of 
the surface abundances is affected by rotation at four dif- 
ferent metallicities. Among the most striking effects one 
can note: 

— The eWNE phase is more extended at high metallicity 
than at low metallicity (see also Fig. inj. This is the 
result of the higher mass loss rates experienced at high 
metallicity. The removal of the outer layers by stellar 
winds is more rapid, therefore when the star enters the 
eWNE phase it does so at an earlier stage of the core 
He-burning phase, when the core is less massive. Thus 
the He-rich envelope extends over a larger fraction of 
the star. 

— For the same physical reason as above, the star en- 
ters the WC phase at an earlier stage of the core He- 
burning phase at high metallicity. Therefore at high 
metallicity one expects that the surface abundances 
will be characterized by higher He abundance and 
larger C/0 ratios than at low metallicity. This has 
interesting consequences for the WC populations ex- 
pected at various Z (see below). 

— The lower the metallicity, the less mass is lost by stel- 
lar winds and thus the less evolved is the stage at 
which the evolution ends. The 40 Mq stellar model 
at Z=0.004 never enters the WC phase. This model 
will explode as a core collapse type II supernova. 

Fig.Elpresents the evolution of the surface H-content 
as a function of the luminosity. The tracks go downwards 
in this diagram. Firstly there is an initial brightening with- 
out surface H-depletion, then due to the concomitant ef- 
fects of mixing and mass loss, the surface abundance of 
hydrogen decreases. The most striking effect of the metal- 
licity is seen for the 120 Mq stellar models. At higher 
metallicity, one notes that the high mass loss prevents the 
track to go through the zone characterized by high lumi- 
nosity and small surface hydrogen abundance. Only the 
models at low metallicity explore this domain of the Xg 
versus Log L/ Lq plane. 

The evolution of the (C-|-0)/He ratio as a function of 
luminosity is shown in Fig. 1161 which is the key diagram 
for WC stars as shown by Smith & Maeder H1991|l . The 
track for a given mass goes up to the left. Indeed as a func- 
tion of time, surface abundances characteristic of a more 
advanced He-burning stage appear at the surface resulting 
in an increase of the (C-|-0)/He ratio. At the same time 
mass loss reduces the actual mass of the star and therefore 
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Fig. 15. Evolutionary tracks in the Xg versus log L/Lq 
plane, where Xg is the hydrogen mass fraction at the sur- 
face. The initial masses are indicated. Long-short dashed 
curves show the evolution oi Z = 0.004 models, dashed- 
dotted curves, short-dashed curves and continuous lines 
show the evolutions for Z = 0.008, 0.020 and 0.040 respec- 
tively. 

its luminosity. Since in rotating models there is a very pro- 
gressive change of the surface abundance, all the models 
would seem to enter the WC phase with more or less the 
same (C-|-0)/He ratio. But this would be misleading since 
this transition phase is very short. Thus we decided here 
to plot the evolution of the surface ratios corresponding 
to the last 90% of the WC lifetime, so that the figure can 
give a better idea of the types of WC stars formed at the 
different metallicities. 

One sees that at high metallicity, mainly late type WC 
stars are present, while at low metallicity, the late type 
WC stars are absent and earlier types are present. This is 
quite consistent with the observational fact that late type 
WC stars are only found in high metallicity regions (Smith 
& Maeder [TMT|l . Fig.lHalso shows that a given WC sub- 
type is reached at lower luminosity for a higher metallic- 
ity, a fact which is also consistent with the observed trend 
(Crowther et al. l2002fl . Such a behaviour results from the 
smaller mass at the entry in the WC phase for models 
at higher metallicity. The smaller mass in its turn is due 
mainly to the much longer WN phase experienced by the 
high metallicity stars, during which intense stellar winds 
peel off the star. 

8. Conclusion 

Rotation gives much better results than non-rotating 
models for the following observed features: 



Fig. 16. Evolution of the ratios (C-|-0)/He as a function 
of the luminosity at the surface of 60 and 120 Af© ro- 
tating models for various initial metallicities (see text). 
Long-short dashed curves show the evolution of Z = 0.004 
models, dashed-dotted curves, short-dashed curves and 
continuous lines show the evolutions for Z = 0.008, 0.020 
and 0.040 respectively. The correspondence between the 
(C-|-0)/He ratios and the different WC subtypes as given 
by Smith & Maeder H1991|l is indicated on the right of the 
figure. 

— the observed number ratio of WR to 0-type stars for 
metallicities between 0.004 and 0.040, 

— the observed ratio of WN to WC stars for metallicities 
lower than solar, 

— the observed fraction of WR stars in the transition 
WN/WC phase 

— and the observed ratio of type Ib/Ic to type II super- 
novae at different metallicities. 

Only in the case of the WC/WN ratio observed at high 
metallicity there may be a difficulty, although complete- 
ness problems may be part of the problem. 

Interestingly, the features which were already well re- 
produced by non-rotating models, such as the abundance 
of ^^Ne at the surface of WC stars or the fact that late 
type WC stars are only observed at high metallicity, are 
also well accounted for by the present rotating models. 
This work further confirms that stellar rotation is an es- 
sential ingredient of massive star evolution, both for the 
MS phase and for the advanced stages like the WR stages 
and the Supernovae. We may therefore also anticipate that 
the nature of the remnants will be different depending on 
rotation and metallicity. In forthcoming papers we shall 
use the present rotating models to estimate the contribu- 
tion of the WR stars to the synthesis of Al in the Galaxy, 
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and also to address the question of the rotation rates of 
pulsars and the progenitors of the coUapsars suggested by 
Woosley l)1993f) to be possible progenitors of long gamma 
ray bursts. 
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